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ABSTRACT 

We present Spitzer/IRAC and MIPS 24 /iin archival observations for 28 brown dwarfs in the Upper 
Scorpius (UppSco) region. We find a disk fraction of 10.7%^g'g%. One object shows a small excess 
at 24 /im but none at shorter wavelengths, and may be a candidate transition disk. Three objects 
show emission in the 10/im silicate feature and we present compositional fits for these sources. Flat 
structures are observed for all disk sources in UppSco. Also presented are the MIPS/70/im observations 
for the TW Hydrae Association brown dwarf 2MASS J1139511-315921. We discuss the structure and 
chemistry of brown dwarf disks at ages of ~5-10 Myr, and consider the possible effects of the brown 
dwarf densities in these clusters on the disk lifetimes. 

Subject headings: circumstellar matter - stars: low-mass, brown dwarfs - open clusters and associ- 
ations: individuals (Upper Scorpius, TW Hydrae Association) - stars: individual 
(2MASS J1139511-315921) 



1. INTRODUCTION 

The Upper Scorpius (UppSco) OB association is part 
of the Scorpius Centaurus complex (de Geus et al. 1989), 
located at a distance of 145±2 pc from the Sun (de Brui- 
jne et al. 1997). The age of USco is -5 Myr, with lit- 
tle scatter of ±3 Myr (e.g., Preibisch & Zinnecker 1999; 
Slesnick et al. 2008). The region is relatively free of 
extinction (Av < 2 mag) and star formation has already 
ended (Walter et al. 1994). The association has been tar- 
geted at multiple wavelengths over the past decade, start- 
ing with the X-rays (e.g. Preibisch et al. 1998) and with 
Hipparcos (de Zeeuw et al. 1999), as well as in the optical 
(Preibisch et al. 2001; Preibisch et al. 2002; Ardila et al. 
2000; Martin et al. 2004; Slesnick et al. 2006; Slesnick 
et al. 2008) and the near-infrared (e.g., Lodieu et al. 
2006; Lodieu et al. 2007). Several tens of brown dwarfs 
with spectral types later than M6 have been confirmed 
as genuine spectroscopic substellar members in UppSco 
(Ardila et al. 2000; Martin et al. 2004; Slesnick et al. 
2006; Slesnick et al. 2008; Lodieu et al. 2006; Lodieu et 
al. 2008). 

Scholz et al. (2007; hereafter S07) conducted a Spitzer 
mid-infrared survey of 35 confirmed sub-stellar members 
of UppSco, selected from the works of Ardila et al. (2000) 
and Martin et al. (2004). Among these 35 sources, 13 
were found to show mid-infrared excesses explained by 
the presence of a circumstellar disk, resulting in a disk 
fraction of 37%±9%. Based on their ^8-24/xm observa- 
tions, S07 noted that objects that do not show any excess 
emission in the mid-infrared are either diskless, or have 
a >5 AU inner hole in the disk. Among the 13 disk- 
bearing objects, only 3 sources were found to show any 
detectable emission in the 10 /im silicate feature. The 
low extent of flaring observed for these 13 sources, along 
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with the flat silicate spectra observed for nearly 75% of 
the objects thus indicated that the UppSco brown dwarfs 
are in a more advanced stage of disk evolution, compared 
to the younger Taurus objects (S07). 

In this work, we consider the optical (RI) survey com- 
plemented by 2MASS infrared photometry, conducted by 
Slesnick et al. (2006; hereafter S06), over the full associ- 
ation down to R ~ 20 mag (corresponding to masses of 
~0.02 Mq). This is the most complete survey available 
to date of the entire UppSco region down to this depth. 
A large number of photometric candidates were selected, 
and a total of 43 were confirmed as spectroscopic mem- 
bers, based on strong Ha emission and weak gravity- 
sensitive features. Here, we present Spitzer observations 
for 28 UppSco brown dwarfs from the S06 survey, along 
with MIPS/70/mi observations for the TW Hydrae As- 
sociation (TWA) brown dwarf 2MASS J1139511-315921 
(2M1139). We compare the structure and chemistry of 
the disks among brown dwarfs at ages of 5-10 Myr, con- 
sidered to be the typical timescale for disk dissipation 
and planet formation. 

2. OBSERVATIONS 

We have analyzed the archival Spitzer/IRAC and 
MIPS 24/mi data (PID 20103; PI: L. Hillenbrand) for 28 
UppSco brown dwarfs (spectral type >M6) from the S06 
survey. S06 have identified a total of 30 brown dwarfs in 
their survey, though archival observations are available 
for 28 of these sources. Aperture photometry was per- 
formed on the post-BCD pipeline data (pipeline versions 
S18.7.0 for IRAC and S16.1.0 for MIPS), using the task 
PHOT under the IRAF package APPHOT. In the IRAC 
bands, the aperture radii ranged between 6 and 12 pix- 
els, depending on the brightness of the target. The inner 
radius of the background annulus ranged between 3-8 
pixels plus the aperture radius, while the width of the 
annulus ranged between 3 and 7 pixels. The inner ra- 
dius and width for the background annulus was changed 
in order to avoid contamination from a nearby source. 
The zero point fluxes of 280.9, 179.7, 115.0 and 64.1 Jy 
were used for the IRAC channels 1 through 4, respec- 
tively. The aperture corrections as noted in the IRAC 
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Data Handbook were applied. The errors in magnitudes 
are between ~0.02 and O.lmag in the IRAC bands. 

In the MIPS/24/mi band, only four out of the 28 Upp- 
Sco targets were clearly detected (detection level of >2a). 
Aperture photometry was conducted using aperture radii 
of 10 pixels for the brighter objects (SCH162630 and 
SCH163247) and 5 pixels for the fainter ones (SCH160930 
and SCH161038), with background annulus of 10-15 pix- 
els and 5-10 pixels for the brighter and fainter sources, 
respectively. A zero point flux of 7.14 Jy was used, and 
an aperture correction of 2.05 and 2.8 was applied to the 
bright and faint sources, respectively. There are three 
sources, SCH162135, SCH162351 and SCH162354, for 
which only MIPS/24/im observations are available. All 
three of these objects were undetected at 24 /im. For 
four sources, low-resolution Spitzer/IRS spectra are also 
available, and these were extracted and calibrated us- 
ing the Spitzer IRS Custom Extraction (SPICE) software 
provided by the SSC. Further details on the processing 
of the IRS spectra are provided in Riaz & Gizis (2008). 
For 2M1139, MTPS/70/im observations were obtained as 
part of the General Observer program (PID 40922; PI: 
Riaz). We had requested an exposure time of 10s and 
20 template cycles. This brown dwarf is undetected at 
70/xm, and we have estimated a 2a upper limit (Section 
%TTj) . The IRAC and MIPS photometry for all sources 
is listed in Table 1. The calibration uncertainty in the 
IRAC and MIPS bands is 2% and 10%, respectively. 

3. RESULTS 
3.1. Spectral Energy Distributions (SEDs) 

Out of the 28 sources in UppSco, the SEDs for 21 ob- 
jects are found to be purely photospheric, i.e. the ob- 
served emission in the IRAC and MIPS/24 /im bands 
coincides with the photospheric model fits within the ob- 
servational errors. To fit the stellar photosphere, we have 
considered the NextGen models (Hauschildt et al. 1999) 
for solar metallicity and log g — 3.5. The effective tem- 
perature was obtained from the spectral type (SpT) using 
the relation for young sources in Luhman et al. (2003). 
The SpT for these brown dwarfs are between M6 and 
M8 (S06). We therefore considered T eff between 3000K 
and 2700K. The uncertainty on the temperature is of the 
order of ~100 K. We have normalized the models to the 
J-band, as it best represents the photospheric emission 
from the central source (e.g., Cieza et al. 2005). 

There are three other sources as noted in Section 
for which IRAC observations are not available, but these 
were undetected at 24 /im. Fig.[T]a, shows the IRAC color- 
color diagram (ccd) for the 25 UppSco sources with IRAC 
observations available. We have overplotted the bound- 
aries considered by Hartmann et al. (2005) to separate 
the Class II and Class III systems. There are 22 objects 
that lie in the lower left region characterized by nearly 
zero excess in both colors, and thus can be classified as 
Class III systems. Hartmann et al. have also found the 
boundary at [5.8]-[8]~0.4 to accurately separate the clas- 
sical T Tauri stars (CTTS) from the weak-line T Tauri 
stars (WTTS). The 22 photospheric objects thus classify 
as WTTS, as is also evident from Fig. [TJd where all of 
these objects show weak emission in the Ha line, and 

4 Chapter 5, Table 5.7 



lie below the empirical accretor/non-accretor boundary 
(S06). The empirical accretor/non-accretor boundary 
has been obtained from Barrado y Navascues & Martin 
(2003). The SpTs and Ha equivalent width measure- 
ments are from S06. 

In the IRAC ccd, three sources, SCH161038, 
SCH162630 and SCH163247, can be classified as Class 
II/CTTS systems. This is consistent with their Ha emis- 
sion being higher than the empirical boundary in Fig.[T]3, 
at a SpT of M6. The SEDs for these brown dwarfs (Fig.Q] 
c-e) show excess emission in at least two of the IRAC 
bands, and at 24/im. One object, SCH160930, shows a 
small excess at 24/im, with a F b s /F p hot of 2.88 (Fig.[TJ). 
The emission at shorter wavelengths, however, is photo- 
spheric. This therefore may be a candidate transition 
disk, i.e. a system with an inner hole in the disk. The 
inner disk region may be either devoid of any disk ma- 
terial or may contain optically thin dust. SCH160930 
is classified as a Class III/ WTTS source, based on its 
IRAC colors and weak Ha emission. Longer wavelength 
observations are required to confirm the presence of a 
transition disk around this brown dwarf. 

We have also considered the classification scheme by 
Lada et al. (2006), which is based on the 3.6-8 fxm slope, 
a, of the IRAC SED and classifies objects into optically 
thick disks, anemic (transition) disks or photospheric 
sources. Sources with optically thick disks are charac- 
terized by a > -1.80, diskless sources by a < -2.56, while 
sources with -2.56 < a < -1.80 are considered anemic 
disks, i.e. objects with optically thin inner disks or disks 
with inner holes. For the three confirmed disk-bearing 
objects SCH161038, SCH162630 and SCH163247, we find 
a of -1.52, -1.31 and -1.67, respectively, thus classifying 
all three to possess optically thick disks. The values of a 
for SCH161038 and SCH163247 are though close to the 
optically thick/anemic disk boundary, which can be ex- 
plained by the weak or no excess emission observed for 
these sources shortward of ^5 /xm. Based on their ^8-24 
/zm observations, S07 had noted that objects that do not 
show any excess emission in the mid-infrared are either 
diskless, or have a >5 AU inner hole in the disk. They 
however could not constrain inner holes of sizes <1 AU 
due to the absence of observations in the ^3-8/im wave- 
length range. For the three sources studied here that 
show excess emission in the IRAC bands, we estimate 
inner disk holes of >0.1 AU (Section §3.2|) . For the can- 
didate transition disk SCH160930, the inner disk radius 
could be as large as ~1 AU. We find an a of -2.6 for 
this source, consistent with it being photospheric in the 
IRAC bands. 

For the case of the TWA brown dwarf 2M1139, Riaz 
et al. (2006) had reported an excess emission at 24/zm 
with a F ba/Fphot of ~3, while Riaz & Gizis (2008) had 
noted an increase in the flux densities in the IRS spec- 
trum for this object, longward of 24/im. The IRS spec- 
trum though was at a very low S/N (<5). Based on the 
same observations, Morrow et al. (2008) had reported 
this object to be photospheric upto wavelengths of ~20 
/im. We have estimated a 2cr upper limit at 70 /im for 
2M1139 by measuring the flux in the regions that lie close 
to the position of the source, but which do not contain 
any other sources. The measured flux thus represents 
the noise level in that region (e.g., S07). Though the 2a 
upper limit for 2M1139 is well in excess of the predicted 
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photospheric emission (F b s /F p h t ~ 70) (Fig. [Tg), this 
brown dwarf is undetected at 70 /*m, and therefore our 
observations are inconclusive regarding the presence or 
absence of a disk around this source. 

3.2. Disk Modeling 

We have used the 2-D radiative transfer code by Whit- 
ney et al. (2003) to model the disks around SCH161038, 
SCH162630 and SCH163247. The circumstellar geome- 
try consists of a rotationally flattened infalling envelope, 
bipolar cavities, and a flared accretion disk in hydrostatic 
equilibrium. The disk density is proportional to vj~ a , 
where w is the radial coordinate in the disk midplane, 
and a is the radial density exponent. The disk scale 
height increases with radius, h = ho(w/R^ , where ho 
is the scale height at i?» and (3 is the flaring power. 
For the stellar parameters, we have used a T e //=3000K, 
M*=0.06 M Q and i?*=0.6 R Q . The effective tempera- 
ture was obtained from the T e f y-SpT relation in Luhman 
et al. (2003), while the mass and radius were obtained 
by considering the 5 Myr isochrone from Baraffe et al. 
(2003) for a T eff of 3000K. The NextGen (Hauschildt et 
al. 1999) atmosphere file for a T e ff of 3000K, and log 
g = 3.5 was used to fit the atmosphere spectrum of the 
central sub-stellar source. A distance of 145±2 pc (de 
Bruijne et al. 1997) was used to scale the output fluxes 
from the models to the luminosity and distance of the 
three brown dwarfs. 

Figs, lc-e show the model fits obtained for SCH161038, 
SCH162630 and SCH163247. Also indicated are the sep- 
arate contributions from the disk, the stellar photosphere 
and the scattered flux. Detailed discussions on the fit- 
ting procedure and the variations in the model SEDs with 
the different disk parameters are provided in Whitney et 
al. (2003) and Riaz & Gizis (2007). Here we provide a 
brief description of the best model-fits (based on the low- 
est reduced-x 2 value) obtained for these brown dwarfs. 
The models used here consider three grain models: large 
grains with a size distribution that decays exponentially 
for sizes larger than 50 /zm up to 1 mm, grains of sizes 
with a max ~ 1 /im, and ISM-like grains with a max ~ 0.25 
fim. We have varied these grain models in the disk mid- 
plane and the upper atmosphere. The object SCH162630 
shows more emission in the 10 /jm silicate feature com- 
pared to the other two sources (Fig.[2f). For this source, 
sub-micron sized grains were placed in both the midplane 
and the disk surface layers. This produces a model SED 
with a more peaked emission at 10/im and a steeper slope 
at far- infrared/sub- mm wavelengths (e.g., Riaz & Gizis 
2007). SCH163247 shows weaker emission while the sili- 
cate feature is nearly flat for SCH161038 (Fig.[2ji and g). 
For both of these, the large grain model with a max of 1 
mm was considered in the disk midplane, while a max ~ 1 
/im grains were placed in the upper layers. This results 
in a weaker silicate emission feature and a flatter sub- 
mm slope in the model SED. For SCH161038, using the 
large-grain model in the disk atmosphere, which results 
in a flat silicate feature, also provides a good fit to the 
observed SED. Considering that we do not have any far- 
infrared/sub- mm observations for these objects, the 10 
/tm feature provides the only way to constrain the differ- 
ent grain sizes considered by these disk models. A more 
detailed compositional analysis of the silicate features is 
provided in Section i )3.3l 



As mentioned, the disk scale height varies as h — 
/io(w/-R*) /3 , where ho is the scale height at i?» and (3 is 
the flaring power. We varied both ho and p to determine 
the amount of flaring in these disks. The variations in 
the model SEDs as (3 is lowered are more evident for A > 
10 /xm, and the 24 /jm observation provides a good con- 
straint to the flaring power. The model SED is highly 
flared for /3=1.25, which is considered a typical value 
for models of T Tauri disks in hydrostatic equilibrium 
(e.g., Walker et al. 2004), while the structure is flat for 
/3=1.0 or smaller values (e.g., Riaz & Gizis 2007). The 
24/im emission for SCH163247 is higher than the other 
two objects, and a slightly higher flaring power of 1.17 
provides a good fit to the observed SED. In comparison, 
SCH161038 and SCH162630 show flatter structures, with 
(3 of 1.1 and 1.0 providing good fits, respectively. The 
scale height at i?* was set to a small value of 0.01 for 
SCH161038 and SCH163247, while a higher value of 2.8 
was required to fit the stronger emission observed near 
10 /im and at shorter wavelengths for SCH162630. None 
of the sources therefore show evidence for strong flaring 
in the disk, as indicated by the low values of j3 between 
1.0 and 1.17. This is consistent with the findings of S07 
who reported a low degree of flaring in the 13 UppSco 
brown dwarf disks studied in their survey. 

We have explored a range of inclination angles to the 
line of sight. Due to binning of photons in the mod- 
els, there are a total of 10 viewing angles, with face-on 
covering 0-18° inclinations. The mid- and far-IR fluxes 
increase with decreasing inclinations, as the emission at 
these wavelengths is from an optically thick disk, while 
the optically thin millimeter fluxes are independent of the 
inclination (e.g., Walker et al. 2004). An inclination that 
is more face-on (<40°) provides a good fit to all three ob- 
jects, with the bin of 18°-32° providing the best-fit. The 
disk outer radius was fixed at 100 AU, since variations in 
this parameters are more evident at far- infrared/sub- mm 
wavelengths (e.g., Riaz & Gizis 2007). As mentioned in 
Section SCH161038 and SCH163247 show photo- 

spheric emission upto wavelengths of ~5/zm, indicating 
inner holes in the disks. For these two sources, an inner 
disk radius of 3R su b was considered, where R su b is the 
dust sublimation radius (lR su b = 3.5i?»). A 3R su b radius 
thus corresponds to an inner hole of ~0.2 AU in the disk. 
For the third source SCH162630, we considered an inner 
disk radius of 1 R su b, equal to about 0.06 AU. Thus for 
the three sources studied here that show excess emission 
in the IRAC bands, we estimate inner disk holes of >0.1 
AU. The disk mass for modeling was set to 1E-4M Q ; this 
is the typical disk mass obtained by Scholz et al. (2006) 
through modeling the 1.3mm observations for a sample of 
Taurus brown dwarfs. For the case of SCH162630, how- 
ever, a steeper slope is observed between the red edge of 
the IRS spectrum and the 24/im point (Fig.[TJi). Reduc- 
ing the disk mass to 3E-7M Q provides a better fit to the 
observed SED for this source. However in the absence of 
millimeter data, it is impossible to constrain the model 
SEDs and obtain a good estimate on the disk mass. 

There are degeneracies in the model fits presented here, 
considering that we have only the IRAC and MIPS/24 
/im data points to model these disks. The uncertain- 
ties in the stellar parameters could also result in changes 
in the disk parameters. There are mainly six free pa- 
rameters related to the disk emission in these models, 
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that can be varied to obtain a good fit. We had fixed 
the disk mass and the outer disk radius given the ab- 
sence of longer wavelength data. There are no mea- 
surements available for the disk mass accretion rates for 
these sources. We had fixed this parameter at a value of 
10 MqJ/T" , a typical value observed among accret- 
ing T Tauri stars (e.g., Muzerolle et al. 2003), consid- 
ering that these sources are more likely to be accreting 
(Fig. []})). However, it is possible to constrain to some 
extent the amount of flaring in the disk, the inclination 
angle as well as the inner disk radius, given the data at 
hand. A wide spread is observed at 24 fim in the model 
SEDs for the different flaring powers, and variations in 
the model SEDs with the disk inclination angles are more 
evident at wavelengths between ~8 and 60 /zm, while 
the optically thin millimeter fluxes are independent of 
the inclination (e.g., Riaz & Gizis 2007). As mentioned 
above, inclinations less than —40° provide better fits to 
all three sources. We can estimate a range in flaring 
power of 1.0-1.1 for SCH161038 and SCH162630, while 
a slightly higher range of 1.15-1.18 provides a good fit 
to SCH163247. The inner disk radii can be estimated 
between —0.1-0.3 AU for all three sources, although a 
radius of 7R su b (— 24i?»«0.6AU) results in higher fluxes 
near the 10 fim silicate band and at longer wavelengths, 
while the emission around 3 fim decreases. This large a 
radius is not a good fit to the observed SEDs. It is thus 
possible to obtain some estimate on the extent of inner 
disk clearing in the disk with the observations available. 

3.3. The lOfim Silicate Emission Feature 

In Riaz (2009; hereafter R09), we had presented 
compositional fits to the 10/im silicate emission features 
for three UppSco brown dwarfs: UScoCTIO J160026.6- 
205632.0 (uscoll2), UScoCTIO J160958. 5-234518. 6 
(usdl60958) and UScoCTIO J161939.8-214535.1 
(usdl61939). The coordinates for these sources are from 
Ardila et al. (2000). To obtain a compositional fit, five 
dust species were considered: amorphous olivine and 
pyroxene, crystalline enstatite, crystalline forsterite, and 
silica. For the amorphous species, two grain sizes of 0.1 
and 2 fim. were considered, while only the sub-micron 
size was considered for the crystalline silicates. A model 
spectrum was then constructed by calculating the sum 
of the emission from a featureless power-law continuum, 
and the emission from the different dust species (at the 
two different grain sizes). To fit the model spectrum, the 
^-minimization method as outined in van Boekel et al. 
(2005) was used. A detailed discussion on the method 
and analysis is provided in R09. The uncertainties 
for the mass fractions were obtained by taking into 
consideration both the noise in the spectrum, and the 
errors for the observed flux densities. In addition to a 
model that included all of the five dust species, R09 had 
constructed two more models: model '(a)' that did not 
consider any crystalline silicates, and model '(b)' that 
included large amorphous grains only. The two addi- 
tional models respectively highlighted the importance 
of the crystallization and the grain growth processes 
in the disk. For the sources uscoll2 and usdl60958, 
the silicate emission features are very weak and show 
a low contrast above the continuum (Fig. Hk.-b). R09 
had found all three models to provide a good fit to the 
observed spectra for these sources, with a difference 



in the reduced-x 2 value between the best-fit and these 
additional fits of <0.1. This thus indicated the large 
uncertainties in their model-fits and the derived mass 
fractions. Among the S06 brown dwarfs, there are 
two similar cases of SCH161511 and SCH161038 that 
show weak, low-contrast features. The flat feature 
observed for SCH161511 is consistent with it being 
a photospheric source. Fig. [2] shows the normalized 
continuum-subtracted spectra for these four objects 
with weak features. In order to normalize the silicate 
feature to the continuum, we first approximated the 
underlying continuum by connecting the two end points 
of the spectrum, i.e. the observed fluxes at 8 and 13 /im 
(averaging over 0.2 /im at the end points; Kessler-Silacci 
et al. 2005). The normalized silicate spectrum is then 
defined in units of [{F v — F c )/F c ], where F v and F c 
are the observed and the continuum normalized fluxes, 
respectively. 

Figure [3] shows the compositional fits obtained for 
SCH162630, usdl61939 and SCH163247. Also shown are 
the additional fits for models (a) and (b). The spec- 
trum for SCH162630 is at a higher S/N of -25, and 
the best-fit indicates negligible grain growth in the disk 
but a high crystalline mass fraction (—46%; Table 2). 
The fits obtained for models (a) and (b) are both of a 
poor quality for SCH162630 (Fig. O, with the differ- 
ence in the reduced-x 2 value between the best-fit and 
these additional fits of >0.5. This object has a more 
flat-peaked spectrum compared to usdl61939, a shape 
which is indicative of a higher degree of dust process- 
ing in the disk (e.g., Sargent et al. 2009; R09). The 
spectra for usdl61939 and SCH163247 are noisy (S/N ~ 
8). While the fits for model (b) using large amorphous 
grains are poor for both objects, the fit for model (a) 
results in a difference in the reduced-% 2 value of only 
0.04, indicating crystalline silicates to be unimportant in 
these disks. However, SCH163247 shows a slight peak 
near 11.3 /ira which is indicative of crystalline forsterite, 
and some fraction of this dust species is required in the 
model to better fit the enhanced emission observed at 
this wavelength. The mass fractions and the reduced-x 2 
values for the best-fits are listed in Table 2. 

Figure 2^ compares the 'shape' and 'strength' in the 
10/im silicate emission feature for the UppSco, Taurus 
(R09) and Cha I (Apai et al. 2005) brown dwarfs with 
higher mass T Tauri stars (Pascucci et al. 2009) and 
Herbig Ae/Be stars (van Boekel et al. 2005). The fea- 
ture strength can be estimated from the peak normalized 
flux, Fp ea k, above the continuum, and the shape can be 
estimated by the ratio of the normalized fluxes at 11.3 
and 9.8/im (e.g., Bouwman et al. 2001). The shapes 
and strengths for the UppSco brown dwarfs are compa- 
rable to the Taurus objects, that cluster around F pea k ~ 
1.3 and -F11.3/-F9.8 ~ 1-0, with the features being weaker 
in these objects compared to most other T Tauri and 
Herbig Ae/Be stars. The brown dwarf usdl61939 shows 
a more narrow, peaked feature than the rest, which is 
consistent with a high 74% small amorphous grain mass 
fraction found for this object. We do not observe in Up- 
pSco or Taurus any of the extreme cases as in Cha I that 
show an increase in the feature strength with increasing 
crystallinity. This 'reversal' in the observed shape vs. 
strength trend for the Cha I brown dwarfs was explained 
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by Apai et al. (2005) due to the presence of highly pro- 
cessed dust grains in the disk. Considering that none 
of the Taurus objects show such a behaviour, and the 
UppSco objects mostly have flat features, the reversal 
seems to be occurring at some intermediate age of ^2-3 
Myr. Fig. HJd compares the disk structures for the three 
UppSco brown dwarfs with the ones in the Taurus and 
Cha I regions. Here, we have plotted the flux densities 
at 8 and 13/mi, normalized to the flux at 8/im. The nor- 
malized flux density at 13/im will be higher for a flared 
disk, and will decrease as the disk geometry gets flatter. 
The processing of dust into crystalline silicates as well as 
growth to larger sizes both effect the vertical structure 
of the disk. While a few sources in Taurus and Cha I 
show flared geometries, all three UppSco brown dwarfs 
show flattened disk structures, comparable to the Taurus 
objects MHO 5 and GM Tau that are dominated by pro- 
cessed dust (~20% crystalline mass fractions; R09). The 
brown dwarf SCH162630 with a 46% crystalline mass 
fraction shows a slightly more flattened disk structure 
than the other two objects, though any dependence of 
the crystallinity levels on the extent of dust settling in 
the disk for the UppSco brown dwarfs is difficult to probe 
for such a small sample. R09 had noted a lack of any 
strong correlation between the two processes for the Tau- 
rus brown dwarfs. 

Figs.[4j;-d compare the crystalline and large-grain mass 
fractions among the UppSco, Taurus and Cha I brown 
dwarfs. In the small SpT range probed here, we do 
not find any particular trends for the crystalline or the 
large-grain mass fractions. The crystalline mass frac- 
tions for the three UppSco sources range between ~14% 
and 46%, with a median value of 26%. For 15 brown 
dwarfs in Taurus that show emission in the 10 /im sili- 
cate feature, R09 had reported a 20% median crystalline 
mass fraction, comparable to the UppSco sources. How- 
ever the very small number statistics of 3 for the UppSco 
objects, and the uncertainties in the mass fractions for 
usdl61939 and SCH163247 due to noisy spectra make 
any such comparison difficult. All three UppSco sources, 
on the other hand, have small amorphous grain mass 
fractions of >50% (Table 2), and smaller large-grain mass 
fractions compared to most of the Taurus brown dwarfs 
(Fig. 0Ji). This may be another indication of the ~5 
Myr UppSco objects to be more evolved systems, since 
increasing grain growth and sedimentation could result 
in an increase in the fraction of small amorphous grains 
in the disk surface layers, as the larger grains settle to- 
wards the disk midplane (e.g., Dullemond & Dominik 
2008; Honda et al. 2006). 

We can consider three main categories for the silicate 
emission features: (a) objects with 'flat' spectra around 
10 /im with the emission levels consistent with the un- 
derlying continuum, (b) objects with 'weak' spectra such 
as uscoll2 that may be transitional objects, and (c) ob- 
jects with emission in the silicate feature. Among the 
20 Taurus objects studied by R09, there is a ~25% frac- 
tion of flat features, and a nearly equal ~20% fraction 
of weak spectra (e.g., CFHT-BD-Tau 20). The detection 
rate for emission is then about 55% in Taurus. Com- 
bining this work and the S07 survey, we have a total of 
16 bona fide disks in UppSco. A 62% fraction of these 
show flat spectra (10/16; S07), while about 19% show 
weak spectra (Fig. [2]), and a similar fraction is found 



for the emission spectra. Thus the fraction of flat fea- 
tures has increased by a factor of ^2 over a ~l-5 Myr 
timescale, while the detection rate in emission has de- 
creased by a similar factor. The fraction of weak features 
that may be in transition from emission to flat spec- 
tra seems comparable in this age range. In the TWA 
(age of ~10Myr; Webb et al. 1999), the two confirmed 
disk sources, 2MASSW J1207334-393254 (2M1207) and 
SSSPM J1102-3431 (SSPM1102), show flat silicate spec- 
tra, though slight absorption due to a more edge-on disk 
is observed for 2M1207 (Riaz & Gizis 2008). The overall 
properties of the silicate emission features are thus con- 
sistent with a picture of increasing dust processing with 
age, since such flat spectra are an indication that signif- 
icant grain growth and dust settling has occurred, and 
the disks are in a more advanced stage of evolution (e.g., 
S07; R09). 

3.4. Disk Evolution in Brown Dwarfs 
From this archival survey of 28 UppSco brown dwarfs, 
we estimate a disk fraction of 10J%~^' 7 Jy (3/28), where 
the uncertainties are the la Gaussian limits. If we con- 
sider the one candidate transition disk, then the fraction 
increases to 14.3%+®3^ (4/28). This fraction is lower 
by a factor of ~3 compared to the 37%±9% (13/35) re- 
ported by S07. Combining the two surveys, we have a 
fraction of 27%+^°^° (17/63) for brown dwarf disks in 
UppSco. In the TWA, out of the four confirmed brown 
dwarf members, two have been confirmed as disk-bearing 
objects (2M1207 and SSPM1102; Sterzik et al. 2004; 
Riaz & Gizis 2008), resulting in a 50%±20% disk frac- 
tion. The TWA disk fraction thus could be somewhere 
between 30% and 70%. This suggests a constant brown 
dwarf disk fraction over a ~5-10 Myr timescale, but could 
also suggest an increase in the fraction by a factor of ^2. 
At younger ages, Luhman et al. (2006) have reported a 
40%±13% brown dwarf disk fraction in Taurus (~lMyr), 
while the fractions are 50% ± 17% and 42% ± 13% in the 
^2-3Myr Cha I and IC 348 regions, respectively (Luh- 
man et al. 2005). Among earlier type stars, Luhman 
et al. (2005) have reported disk fractions of 45% ± 7% 
and 33% ± 4%, respectively, in the Cha I and IC 348 re- 
gions for the M0-M6 members (0.7 M > M > 0.1 M ). 
Carpenter et al. (2006) have reported a 19% ± 5% frac- 
tion for the K0-M5 stars in UppSco, while Low et al. 
(2005) have found a disk fraction of 24%+^ % (6/25) for 
the stellar members in the TWA. If we consider only the 
M0-M6 stars, then the TWA disk fraction drops down to 
16.6%^g 1 4 y o (3/18). A simple comparison then suggests 
that the brown dwarf disk fraction has not decreased by 
an age of ^10 Myr, while it shows a drop by a factor of 
~2 among the higher mass stars. Brown dwarf disks thus 
may survive until ^10 Myr, as earlier noted by e.g., S07, 
Bouy et al. (2007) and Riaz & Gizis (2008). The current 
statistics in the TWA are too limited to put any limits on 
the disk evolutionary timescale for brown dwarfs. Also, 
all such conclusions are affected by the uncertainty on 
the age estimate for any given cluster. 

There is one source SCH160930 that may be a candi- 
date transition disk, resulting in a fraction of 3.6% jl! 1 
of such systems among the UppSco brown dwarfs. This 
is consistent within the uncertainties with the 8% upper 
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limit estimated by S07 for brown dwarf transition disks in 
UppSco. SCH160930 shows photospheric emission short- 
ward of 24/im, indicating an inner disk hole as large as 
~1 AU. Such large inner holes have been reported earlier 
for brown dwarfs in the IC 348 cluster (Muzerolle et al. 
2006). The absence of excess emission at shorter wave- 
lengths could be either due to grain removal processes 
such as photoevaporation, or substantial grain growth 
and dust settling in the inner disk regions (e.g., Muze- 
rolle et al. 2006; Cieza et al. 2007). 

The disk fraction obtained here is a factor of ~3 smaller 
than that reported by S07 for brown dwarfs in UppSco. 
In Fig. [5J we have examined the spatial distribution for 
the disk-bearing sources in the S07 sample and in this 
work. While 80% of the S06 sources have A v < 1, the 
source SCH162630 lies in a more extincted region to- 
wards the south-east, and S06 have measured an A v of 
3.02 for this object. SCH161038 that lies in the northern 
region also shows some extinction of A„=1.78 (S06). The 
third disk SCH163247 however lies in a clear region with 
an A v of 0.06. While these three sources lie to the east 
of most disk-bearing objects from the S07 sample, we do 
not see any particular differences in the spatial distribu- 
tion of the disk sources in the two samples, that might 
indicate an age effect. 

We consider the brown dwarf densities in the Taurus, 
UppSco and TWA regions, that might affect the disk 
fractions noted above. According to Briceno et al. (2002) 
and Luhman et al. (2003), the number of brown dwarfs 
in Taurus with SpT between M6 and M9 and A v < 4 
mag is 10 over 8.4 square degrees, implying a density 
of ~ 1 object per square degree. For UppSco, S06 con- 
firmed 43 members over 200 square degrees from a sam- 
ple of 20% of all photometric candidates. Out of those 
43 sources, 30 are brown dwarfs (SpT>M6). Thus the 
total survey should contain about 5x30=150 objects in 
200 square degrees that have masses between 75 and 20 
Jupiter masses, implying a density of 0.75 per square de- 
gree. From the work of Lodieu et al. (2007), there are 23 
brown dwarfs with estimated masses between 75 and 20 
Jupiter masses over 6.5 square degrees, implying a den- 
sity of 3.5 brown dwarfs per squre degree. Their survey 
however is targeted more towards the core of UppSco, 
which could explain the higher densities. The difference 
with Taurus in terms of density is thus not very large 
and of the same order. This is consistent with the com- 
parable disk fractions observed in the two regions (~30% 



in UppSco, ^40% in Taurus). In comparison, the TWA 
is very loose and contains a small number of members 
(e.g., Webb et al. 1999). Following Jayawardhana et al. 
(1999), the number of young stellar objects belonging to 
the TWA in the solar neighbourhood would be of the or- 
der of 0.06 objects per square degree or less. Considering 
a typical value of ~5 for the ratio of stars to sub-stellar 
objects in young nearby clusters (Andersen et al. 2008), 
we can estimate a brown dwarf density in the TWA that 
is at least ^100 times less than that in Taurus or Up- 
pSco. The TWA however is a factor of ^3 closer than 
Taurus or UppSco, and therefore even with the same ob- 
ject density per cubic parsec, we would expect a drop 
by a factor of ~10 in the density per square degree. Al- 
though it can be argued that the larger age of the TWA 
allows more time for the dispersion of the young sources, 
the much lower densities in this region could explain the 
longer disk dissipation timcscalcs, and disk lifetimes may 
be prolonged in a less dense environment. 

4. SUMMARY 

We report a disk fraction of 10. 7%^' ^ from a Spitzer 
archival survey of 28 UppSco brown dwarfs. One object 
shows a small excess at 24 //m but none at shorter wave- 
lengths, and may be a candidate transition disk. The 
fraction of flat silicate spectra has increased by a factor 
of ~2 over a ~l-5 Myr timescale. While a few Taurus 
sources show flaring in the disk, all disk-bearing brown 
dwarfs in UppSco and the TWA show flat disk structures, 
consistent with a picture of increasing dust processing 
and dust settling with age. We do not find any differ- 
ences in the spatial distributions for the disk-bearing ob- 
jects found in this work, and from the S07 survey, that 
might indicate an age spread in the UppSco. The much 
lower brown dwarf densities in the TWA might be an 
explanation for the longer disk lifetimes observed among 
the sub-stellar objects in this region. 
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IRAC AND MIPS MAGNITUDES 
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TABLE 2 

Model Fit Parameters 



Name 


Fpeak 


F11.3/F9.& 


Amorphous a 
Silicates 


Large b 
Silicates 


Crystalline 
Silicates 


x 2d 


SCH161038 
usdl61939 
SCH162630 
SCH163247 


1.56 
1.38 
1.37 


1.0 
0.68 
0.7 


48.5%^ 

OM -'-13.6% 


0.005%™ 
0.004%±4% 
26.6%t^ 


25.9%+ 6 4 ™ 
46.2%+f™ 

13.6%^;^ 


1.6 
1.1 
1.4 



a Percentage of small amorphous olivine and pyroxene silicates compared to all other silicates. 

Percentage of large amorphous olivine and pyroxene silicates compared to all other silicates. 
c Percentage of crystalline silicates (enstatitc and forstcrite) compared to all other silicates. 
d The reduced-x 2 value of the best-fit. 
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FlG. 1. — (a): IRAC ccd for 25 UppSco sources, (b): Ha emission vs. SpT for the 28 brown dwarfs studied here. Dashed line marks the 
empirical accrctor/non-accretor boundary. Objects with disks are denoted by red circles, while the one candidate transition disk SCH160930 
is denoted by a blue circle, (c)-(e): Best model fits for the three disk-bearing objects, (f)-(g): SEDs for SCH160930 and 2M1139. 
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Fig. 2. — The normalized continuum-subtracted spectra in units of (F v — F c )/F c for the weak/flat features. Dashed line represents the 
continuum. 



SCH162630 SCH162630 SCH162630 




usd161939 usd161939 usdl61939 




Fig. 3. — Model-fit to the 10/xm silicate features. Left panel shows the fits obtained using all five dust species, middle panel shows the fits 
obtained without using any crystalline silicates, right panel shows the fits obtained using only large amorphous olivine and pyroxene grains. 
Colors represent the following: red-observed spectrum; black— model fit, grey— small amorphous olivine; green-small amorphous pyroxene; 
cyan— silica; blue-forsterite; pink— enstatite; black dashed— large amorphous olivine; orange-large amorphous pyroxene. Thin dashed line 
represents the continuum. 
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Fig. 4. — (a) The shape vs. strength of the silicate emission features. Taurus brown dwarfs are indicated by filled black squares, UppSco 
ones by filled red circles. Included for comparison are T Tauri stars from Pascucci et al. (2009) and Herbig Ae/Be stars from van Boekel 
et al. (2005). (b): Disk geometry for the Taurus, UppSco and Cha I brown dwarf disks. The values next to the target names are the 
crystalline mass fractions, (c)-(d): Crystalline and large-grain mass fractions vs. SpT. The values of 5-9 indicate SpTs of M5-M9. 
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Fig. 5. — Spatial distribution for the disk-bearing objects from this work and the S07 survey. The image is a false-color 13°x8° IRAS 
12 fim observation. The color intensity scale in units of MJy sr — 1 is shown at the bottom, and the coordinate grid in degrees (J2000) is 
overlaid. Symbols represent: black circles— S06 brown dwarfs without disks; black crosses— S07 objects without disks; blue squares and blue 
diamond represent the three disk-bearing objects and the one candidate transition disk, respectively, found in this work; green squares 
represent the disk-bearing objects from the S07 survey. 



